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Abstract 

In a recent paper [S.El Boustani, P.R.Buenzli, and Ph.A.Martin, Phys.Rev. E 73, 036113 (2006)], 
about quantum charges in equilibrium with radiation, among other things the asymptotic form of 
the electric-field correlation has been obtained by a microscopic calculation. It has been found that 
this correlation has a long-range algebraic decay of the form 1/r^ (except in the classical limit). 
The macroscopic approach, in the Course of Theoretical Physics of Landau and Lifshitz, gives no 
such decay. In this Brief Report, we revisit and complete the macroscopic approach of Landau and 
Lifshitz, and suggest that, perhaps, the use of a classical electromagnetic field by El Boustani et 
al. was not justified. 
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I. INTRODUCTION 



We consider the equilibrium quantum statistical mechanics of an infinite and homo- 
geneous system of nonrelativistic charged particles coupled to electromagnetic radiation. 
For this system, we investigate the correlation function of the electric field at equal times 
< Ei(ji)Ej{v2) > where -Ej(ri) is the quantum operator for the ith Cartesian component 
{i = 1,2,3) of the electric field at point ri and < . . . > denotes a quantum statistical average 
at temperature T; in the present homogeneous medium, this correlation function depends 
on ri and r2 only through r = ri — r2 = {xi, X2, X3). We want to compute the form of this 
correlation function, when r becomes large compared to the microscopic scale. 

A macroscopic approach to this problem has been made a long time ago by Landau and 
Lifshitz (LL) in their famous Course of Theoretical Physics[l|, Actually, their theory has 
been written for any medium, characterized by a complex frequency-dependent dielectric 
function e{uj). Recently, a microscopic theory has been elaborated by El Boustani, Buenzli, 
and Martin (BBM) . They find an electric- field correlation function in disagreement with 
the one advocated by LL: while BBM find that the correlation function has a long-range 
power-law decay of the form (except in the classical limit), I could not extract this 
algebraic decay from the work of Landau and Lifshitz. The reason for this disagreement is 
an open problem. 

As a first step for clarifying this problem, the present paper revisits and completes the 
macroscopic approach of LL. Section II summarizes this approach, making it more explicit 
about the decay of the electric-field correlation function. In Section III, the formalism is 
applied to the special case of a one-component plasma. Section IV is a Conclusion where 
the above-mentioned discrepancy is discussed. 



II. MACROSCOPIC APPROACH 

A. The results of Landau and Lifshitz 



LL have solved the macroscopic Maxwell equations, in a medium characterized by a 
complex frequency-dependent dielectric function e{u!), in presence of a random field. This 
macroscopic approach is expected to be valid only for distances large compared to the 
microscopic scale. The fiuctuation-dissipation theorem is used for obtaining the electric- 
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field correlation as a Fourier transform with respect to r and the time difference (here this 
time difference is zero): 

r~ duj r d^/c 



< Ei{ri)Ej{r2) >= f tt" / ^ exp(ik ■ r)Sij{k, u) 

J~oo ZtT J 071'^ 



(1) 



where, in terms of the wave vector k and the frequency u (see §76 and 77 of 



0: 



£:.,(k,.) = -4vr;.coth^Im ^^^^J^^/2)-F 



A. . U.k- 

uo^t{uo) 



(2) 



h is Plank's constant divided by 27r, c is the velocity of light in vacuum, T is the temperature 
in units of energy. 

In order to make easier the comparison with the paper by BBM, we shall now separate 
([2]) in its longitudinal and transverse parts. 



B. Longitudinal and transverse correlations 

Introducing the projectors kikj/k"^ and 5ij — [kikj/k"^) on the longitudinal and transverse 
parts of the electric-field correlation functions, we can rewrite ([2]) as 

£,,(k,o;) = 4.(k,.;)+4.(k,.;), (3) 

where the longitudinal part is 



4(k,^) = -47r;icoth^Im-^^ (4) 



and the transverse part is 



£.,(k, u.) = -Ann coth [6., - ^ j • (5) 

It should be noted that the cross correlation between the longitudinal and transverse parts 
of the electric field vanishes. Indeed, in Fourier space, in terms of the Fourier trans- 
forms of these fields E^(k, cj) and E*(k, cu), the cross correlation tensor is proportional to 
< E'(k, co')E*(— k, —uj) >, where E' is along k and E* is normal to k. Since the medium is 
isotropic, the correlation tensor is unchanged if E* is replaced by its opposite; therefore, the 
correlation tensor is equal to its opposite, i.e. it vanishes. 

The integral on of (j4]) converges: in particular, near u = 0, the dielectric function e{uj) 
behaves like Ama/oo, where a is the static conductivity [l] (see Section III for the special case 



of a one-component plasma) and (jl]) is finite at u = 0. The Fourier transform of Air/k"^ is 
1/r and the longitudinal part of ([1]) is 
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(6) 



At small /c, in ([5]), the term {u"^ j c^) j [{u"^ / c^)e{ijj) — k"^)] can be replaced by l/eiu). Indeed, 
if 7^ 0, this term can be expanded in powers of k"^ and l/e{ijj) is the leading term. If 
\uj\ < ojq, where ojq is a sufficiently small constant, t{uj) = inicr/u!, /icoth|^ = 2T/u!, and 



'■^0 duj . „ , 2T 
— £^(k,^) = — 



arctan 



[ckf 



(7) 



At small k, the arctan in ([7]) behaves like 7r/2. Finally, for all values of cu, at small /c. 



4(k,^) ~ -47r^coth^ 



Im- 
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The asymptotic behavior of the transverse part of ([T]) is given by the most singular part 
at small of ([8]), which is just opposite to (jl]), i.e. the asymptotic form of < Ei{Yi)Ej{Y2) >t 
is opposite to (Ej). These asymptotic behaviors of the form cancel each other in the total 
correlation function ([T]). This cancellation had been previously noted 0, Q] in the classical 
limit T —>■ 00, but the present macroscopic approach predicts that this cancellation persists 
in the quantum regime at any temperature, contrarily to the prediction of the microscopic 
theory of BBM. 

The integral on in ([6]) is simply related to the second moment of the charge correlation 
function. Indeed the charge density p is given by the Poisson equation divE = Anp. Only 
the longitudinal part of E has a non-vanishing divergence. Therefore, since the Fourier 
transform (jlj) of < Ei{ri)Ej{r2) >i is of the form Akikj/k"^, the Fourier tranform of < 
p(ri)p(r2) > is (An)^'^ Ak"^ , which means that A is — (Svr)^^ times the second moment of this 
charge correlation function. Taking the inverse Fourier transform of Akikj/k"^ gives for the 
longitudinal field correlation 



< EAri)EJr2) >i 



dxidxj r 



_^ f dVr'^ <p(0)p(r') > 

o J 



(9) 
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for r large compared to the microscopic scale, in agreement with a microscopic derivation 

In the classical limit T 00, ^coth ^ ~ 2T /oj and the integral on u in ([6j) (which also 
occurs in the transverse part with the opposite sign) has the simple value — T 6|. This can 



be shownjs] by invoking that e{uj) has no zeros when uj is in the complex upper half-plane. 
The calculation is made in the Appendix. Therefore, in the classical limit, 



and 



< Ei(ri)EJr2) >i~ -T 



< EAv^EAv^) >t~ T 



Q2 



dxidxj r ' 
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dxidxj r ' 



(10) 



(11) 



these classical results can also be obtained from since, in the classical case, the second 
moment in obeys the Stillinger-Lovett sum rulejoj] 

27T 



— I dV r'2 < p(0)p(r') >= T. 
3 J 



(12) 



The asymptotic classical transverse correlation functi on ([TTD is the one of a free field. This is 



in agreement with the Bohr-van Leeuwen theorem 



Id, 
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121 ] which says that, in equilibrium 



classical statistical mechanics, matter and radiation are uncoupled. 



III. ONE-COMPONENT PLASMA 



A. Drude dissipationless dielectric function 



The above formalism simplifies in the special case of a one-component plasma, a system 
of one species of particles of charge e, mass m, and number density n, in a neutralizing 
homogeneous background. We are interested in small wave numbers, when the dissipation 



goes to zero[l3|] (i.e. the static conductivity is infinite), and the dielectric function can be 



taken as the Drude one 



= 1 - , X- y (13) 

uj{uj + irj) 

where cUp = (47me^/m)^/^ is the plasma frequency and the dissipation constant t] is taken 
as infinitesimal. 

Then, using Im[l/e(co')] = —'k{uj/\uj\)uj'^5{uj'^ — cjp) in (jl]), one finds for the longitudinal 
correlation function 

< E.ir,)E,ir,) » = -th ^ (14) 



5 



In a similar calculation, using (fT3|) in ([5]), S^uj"^ — uj"^ — c^k"^) appears, and one finds for 
the transverse correlation function 



< 



/Ait. 1 
|^exp(ik-r)-M^J + cV)V^ 



X coth 47r \6ij - -p^ J • (15) 

For small k, the ^-singularity in (fT5|) comes from the term kikj/k"^ with the replacement of 
(cijp + c^k'^Y^'^ by cijp. Then, one finds for the asymptol 
function just the opposite of the longitudinal one ( fT^ . 



(cijp + c^k'^Y^'^ by cijp. Then, one finds for the asymptotic form of the transverse correlation 



B. Longitudinal and transverse modes 

The above results for the electric-field correlation functions can also be obtained, perhaps 
in a more transparent way, from the modes of vibration of the plasma. In this subsection, 
we shall not use the dielectric function e(c(j), but rather take into account explicitly in the 
equations the charge and electric-current densities. 

For every wave vector k there is a longitudinal mode; its frequency jl^ is uj^ (neglecting 
a term of order fc^, which is consistent with the previous use of an e independent of k). In 
such a mode, the electric field, at position r and time t is of the form 

Ek(r,t) = Re(k/A;)Eoexp(ik-r-icjpt). (16) 

Since this mode is an oscillator, for studying it in quantum mechanics, we can first use 
classical mechanics, and quantize at the end. A collective velocity is given by Newton's 
law: 

= eEk. (17) 

From ( JTTj) one easily finds that the temporal average of the kinetic energy density associated 
to Vk is equal to the temporal average of the energy density of the electric field. Taking into 
account that the temporal average of the squared real electric field is \Eq\^ /2 and equating 
the temporal average of the total energy to the statistical average value for an oscillator 
(including the zero-point energy) gives 

^ = i^coth% (18) 
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where V is the volume of the system. The contribution of this mode (fT6|) to the electric-field 
correlation is Re(|i?oP/2) exp(ik ■ r){kikj/k'^), where l-EoP is given by f|T8|) . Summing on all 
k wave vectors, i.e. computing the integral V J d^A;/(87r^) . . . reproduces fll4l) . 

For every wave vector k, there are also two (there are two possible polarizations) trans- 
verse modes of frequency IJ] cuk = i^l + c^k'^Y^'^. For such a mode, the electric field is of 
the form 

Ek(r,t) = Ree£^oexp(ik ■ r — itUfct), (19) 
where e is a unit vector normal to k, and the magnetic-induction field is 

Bk(r,t) = — kx Ek. (20) 

Again, a collective velocity is given by Newton's law f|T7j) (the magnetic force can be 
neglected for non-relativistic velocities). The temporal average of the electric-field en- 
ergy density is |-E'oP/(16vr). The temporal average of the kinetic energy density is 
(co'p/co'^)|£'oP/(167r), from (IT7|) . The temporal average of the magnetic-induction energy 
density is (c^/c^/co'|)|£'oP/(167r), from fl2U]) . Thus, the temporal average of the total energy 
density again is |i?oP/(87r), and 

coth— -. (21) 



8n V 2 2T 

The contribution of the transverse mode (fT9|) to the electric-field correlation is 
Re ejej(|-EoP/2) exp(ik ■ r) where |-EoP is given by fl2T]) . Adding the contribution of the 
other polarization replaces CiCj by 6ij — [kikj/k'^). Summing on all k wave vectors repro- 
duces f lT^ . 



IV. CONCLUSION 



The macroscopic approach of LL is in agreement with the microscopic calculation of 
BBM, for the longitudinal part of the electric-field correlation (incidentally, this agreement 
is an indication that the macroscopic approach can be correct). The disagreement is about 
the transverse part only, which exactly cancels (jH]) in the macroscopic approach while it 
obeys (ITTl) in the microscopic approach[.3], even in the quantum regime. I am tempted to 
believe that the cancellation predicted by the macroscopic calculation is essentially correct. 
Here are my reasons. 
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BBM[3] point out that the macroscopic theory of LL uses a local dielectric function e{uj) 
rather than a k-dependent one. Indeed it is tempting to use a longitudinal ei{k,uj) in (jl]) 
and a transverse et{k,uj) in (jSj). However, the leading (singular) term of (jlj) or ([2]) at small 
k would be still obtained by taking these fc-dependent dielectric functions at A; = where 
both of them become e{uj). Thus, a singular term 1^ in the small-A; expansion of e\{k,uj) 
and/or et{k,uj) would not change the prefactors of the terms k^kj/k"^ (this cancellation only 
concerns the asymptotic terms of the form however an algebraic decay faster than 
could remain). BBM make another criticism of the macroscopic approach of LL: it neglects 
the magnetic permeability. It would be strange that taking the magnetic permeability into 
account in the macroscopic approach would change the transverse correlation function into 
the one of a free field (however, it might bring small corrections to the results of LL[16]). 

The final remark by BBM might be the key point: in their microscopic approach, they 
use a classical electromagnetic field. That their asymptotic transverse correlation is the one 
of a free field, decoupled from matter, even in the case of quantum particles, might be due to 
this feature of their calculation. In III.B, I have argued that the transverse modes certainly 
feel the presence of matter. Perhaps the argument, in the 3rd paragraph of the Introduction 
of BBM, in favor of using a classical electromagnetic field, has a flaw: The frequency of a 
transverse mode is always larger than the plasma frequency cUp, and the condition f3hck <^ 1 
does not imply that (3?LUjp <^ 1. 

The microscopic calculation of BBM is very elaborate, cleverly using an elegant path- 
integral formalism. It would certainly be very interesting to redo this microscopic calculation 
with a quantized electromagnetic field, if feasible. 

APPENDIX 

This Appendix is based on a private message from B.U.Felderhof. In the classical limit, 
the integral in becomes TI where / is 




(A.l) 



We want to show that / 



-1. 



Since e*{uj) = e(— u;*), / can be rewritten as 




(A.2) 
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Since e^uj) has no zeros in the upper complex half-plane the integral along the real axis 
in flA.2|) can be changed into an integral along the half-circle C at infinity in the upper 
half-plane. Since e is 1 at infinity, this integral is J^du/uj = —itt. Therefore, / = —1. 

The same reasoning applies to the transverse part. Comparing (jlj) and one sees 
that, on the half-circle at infinity, 1/e and {u"^ /c^)/[{uj'^ /c^)e — k"^] have the same limit 1. 
Furthermore, the term kikj/k"^ in (j4]) and the term 5ij — {kikj/k"^) in ([5]) will have opposite 
inverse Fourier transforms for r 7^ 0. Thus, ffTTl) is just opposite to ffTOj) . 
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